Band renormalization of a polymer physisorbed on graphene investigated by 

many-body perturbation theory 
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Many-body perturbation theory at the GoWo level is employed to study the electronic properties 
of poly(para-phenylene) (PPP) on graphene. Analysis of the charge density and the electrostatic 
potential shows that the polymer-surface interaction gives rise to the formation of only weak surface 
dipoles with no charge transfer between the polymer and the surface. In the local-density approx- 
imation (LDA) of density-functional theory, the band structure of the combined system appears 
as a superposition of the eigenstates of its constituents. Consequently, the LDA band gap of PPP 
remains unchanged upon adsorption onto graphene. GqWq calculations, however, renormalize the 
electronic levels of the weakly physisorbed polymer. Thereby, its band gap is considerably reduced 
compared to that of the isolated PPP chain. This effect can be understood in terms of image 
charges induced in the graphene layer, which allows us to explain the quasi-particle gap of PPP ver- 
sus polymer-graphene distance by applying a classical image-potential model. For distances below 
4.5 A, however, deviations from this simple classical model arise which we qualitatively explain by 
taking into account the polarizablity of the adsorbate. For a quantitative description with predic- 
tive power, however, we emphasize the need for an accurate ab-initio description of the electronic 
structure for weakly coupled systems at equilibrium bonding distances. 

PACS numbers: 71.10.-w,73.20.-r,85.65.+h,71.15.Mb 



I. INTRODUCTION 

When a molecule is brought in contact with a sur- 
face, its electronic states are generally renormalized in 
terms of their energy positions and widths, where the 
magnitude of these effects depends on the adsorption 
distance. Molecule-surface interactions may give rise to 
several phenomena including hybridization of molecular 
states with substrate levels, charge transfer between the 
molecule and substrate, and formation of strong short- 
range interface dipoles. All of these effects can be de- 
scribed within a static mean-field approach as entailed 
in density-functional theory by solving single-particle 
Kohn-Sham (KS) equations and employing standard lo- 
cal (LDA) or semi-local (GGA) approximations for the 
exchange-correlation potential. 

An additional effect that can greatly alter the level po- 
sitions of the molecule, but is not accounted for in such 
a mean-field approach, arises due to polarization effects. 
Electron addition (removal) energies of the molecule, i.e., 
the quasi-particle energies of the molecule, inherently in- 
volve a negatively (positively) charged molecule residing 
on top of a polarizable substrate. Thus, electro-static 
Coulomb forces of this added electron or hole polarize 
the underlying surface which, in turn, affects the energy 
position of the added electrons or holes, respectively. In 
particular, the ionization potential and electron affinity 
level, and hence the band gap of the molecule, are renor- 
malized due to the presence of the substrate*^ Even for 
physisorbed and weakly coupled molecules, these surface 



polarization effects can considerably change the size of 
the gap of the molecule compared to its gas phase value. 

The response of the electronic system to the added 
electron/hole is not captured by local or semi- local KS 
functionals. In order to account for this dynamical polar- 
ization effect, we employ many-body perturbation the- 
ory on top of DFT calculations. In the GqWq approx- 
imation, the self-energy £ is given as the product of 
the noninteracting single-particle Green function, Go, 
and the screened Coulomb interaction, Wo, calculated 
within the random-phase approximation^ In this ap- 
proach, correlation effects including the response of the 
electronic system to the added electrons or holes are 
taken into account by the nonlocal and energy-dependent 
self-energy. As already demonstrated for the case of a 
benzene molecule physisorbed on various metallic and 
semi-conducting substrates^ the GoWo band gap of 
the molecule is reduced relative to its corresponding gas- 
phase value by an amount that depends on the polariz- 
ability of the surface. In contrast, DFT calculations em- 
ploying LDA or hybrid exchange-correlation functionals 
like PBEO, render the size of the benzene HOMO-LUMO 
gap independent of the underlying substrate. 

In this paper, we explore how the electronic structure 
of a one-dimensional system - the polymer poly(para- 
phenylene) (PPP) used as active material in blue light- 
emitting diodes^ - is affected by the presence of the two- 
dimensional graphene layer. We thereby extend previ- 
ous ab initio worki^ which only treated small organic 
molecules on various surfaces. We study the renormal- 
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ization of the molecular electronic levels by employing 
G W calculations. We first provide details on our com- 
putational methodology concerning the DFT and Go Wo 
calculations. Then, the electronic band structure of 
graphene is discussed, both in the hexagonal unit cell 
and in a larger rectangular supercell necessary for the cal- 
culations of the adsorbate system. The electronic band 
structure at the LDA and Go Wo level of isolated PPP is 
presented and compared to available literature results. 
As the main results, we discuss the combined system 
(denoted as PPP@gr) showing that the polymer's quasi- 
particle HOMO-LUMO gap is strongly renormalized due 
to the presence of the graphene layer. We model our ah 
initio results by a classical image-potential model in or- 
der to describe the asymptotic behavior of the potential 
felt by an electron outside the substrate. We also extend 
the standard image-charge model as to include the polar- 
izability of PPP leading to a better model-description of 
the Go Wo results at short polymer-graphene separations. 

II. COMPUTATIONAL APPROACH 



the calculation of Wo. In this work, the frequency de- 
pendence of Wo is described by a plasmon-polc model 
(PPM) 3 



In the calculation of the QP energies, there are two 
main technical problems which make the GoW) approach 
computationally costly for large unit cells as necessary 
for studying the adsorption of molecules or polymers on 
surfaces. The first bottleneck arises from the summation 
over unoccupied states which appears in the correlation 
part of the self-energy. To overcome this problem, we 
make use of the recently developed energy effective tech- 
nique (EET)pii In this method, all necessary steps in a 
GoWo calculation can be restricted to occupied states 
only. As demonstrated in the Appendix for graphene 
(Fig. [9]), this approximation preserves the precision of the 
conventional sum-over-empty-states approach but speeds 
up calculations by more than an order of magnitude. 



Self-consistent DFT calculations are performed by us- 
ing ABINITpl which is a plane- wave based code with 
periodic boundary conditions along all three directions. 
We utilize norm-conserving pseudo-potentials generated 
by the Troullier-Martins schemed Exchange-correlation 
effects are treated within the local-density approxima- 
tion. An energy cutoff of 50 Ry is chosen for the elec- 
tronic wave functions. Brillouin zone integrations are car- 
ried out by using Monkhorst-Pack^ meshes of (5x17x1) 
k points for the PPP@graphene supercell along with a 
Mcthfcsscl-Paxton^ smearing with a smearing parameter 
of 0.01 Ry. The so obtained DFT ground-state results 
serve as a starting point for the subsequent Go Wo com- 
putations. 

In the GoW]@LDA method, the QP energies are ob- 
tained from the linearized QP equation: 



T QP 
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Here, i>n DA and e^ DA are DFT-LDA eigenstates and 
eigenvalues, and the rcnormalization factor Z n is given 
by 
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The self-energy, E, is calculated non-self-consistcntly 
from the convolution of the non-interacting single- 
particle Green function G with the screened Coulomb in- 
teraction Wo, i.e., £ = iGoW]. Here, the Green function 
is defined as Gq(z) = (z — H LDA )^ 1 where the subscript 
symbolizes the fact that these quantities are obtained 
from the Kohn-Sham orbitals and energies in a non-self- 
consistent manner. Vertex corrections are neglected both 
in the self-energy and in the polarizability and, hence, in 



The second difficulty originates from the periodic 
boundary conditions in combination with the repeated- 
slab approach to model surfaces. Thereby, one com- 
monly constructs unit cells containing a vacuum layer 
to separate periodic images in order to avoid spurious 
interaction between them. In ground-state calculations, 
this poses no big problem since electronic wavefunctions 
and the (semi) local exchange-correlation potential decay 
rapidly into the vacuum, and a moderate thicknesses 
of the vacuum slab suffices. For GW calculations, the 
main difficulty arises from the nonlocal nature of the 
self-energy, in particular, nonlocal Coulomb matrix ele- 
ments, which lead to long-range image charge effects that 
converge very slowly as a function of the inserted vac- 
uum layer To overcome this obstacle, we make use of 
a truncated Coulomb potential which prevents the inter- 
action between periodic images. We utilize Ismail-Bcigi's 
method^ for the case of sheet-like geometries with one 
confined and two periodic directions. The truncation 
length is chosen to be half the lattice parameter per- 
pendicular to the surface (in z direction). Using this 
method, Go Wo band energies converge fast as a function 
of vacuum layer thickness as can be seen from Fig. 1101 in 
the Appendix. Care must, however, been taken in the 
choice of the fc-grid since the convergence with respect to 
the number of fc-points is somewhat slowed when using a 
cut-off Coulomb potential compared to the convergence 
behavior for the plain Coulomb potential (which will be 
seen from the right panel of Fig. [TU)) . As in the public 
version of ABINIT, the truncation of the Coulomb inter- 
action for slab geometries is not fully implemented, we 
have added one term taking the q — > limit. This way, 
we could assure fast convergence of results with respect 
to fc-grid and vacuum size. 
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A. Graphene g 



Before we presents results for the adsorbed polymer 
on graphene (compare Fig. [5]) , we review the band struc- 
ture of uncovered graphene at the LDA and Go Wo level. 
This is depicted in Fig. [TJ where the LDA (Go Wo) band 
structure is shown as dashed (solid) lines. As already 
noted earlier^ the GoWo corrections slightly enhance 
the overall band widths of tt and a bands; for instance, 
they increase the gap at the T point from 6.6 eV to 7.0 eV. 
Close to the Dirac point (K) , the Fermi velocity resulting 
from the LDA dispersion is about 1.01 x 10 6 ms" 1 while 
the Go Wo value of 1.11 x 10 6 ms _1 is about 10% larger 
and improves the agreement with experiment (1.1 x 10 6 
ms -1 ) J£ These results emphasize the fact that the Go Wo 
quasi-particle band structure not only provides improved 
band gaps but also leads to an improved description of 
graphene in terms of band widths and Fermi velocities. 

In order to study the adsorption of PPP on graphene, 
we need to construct an appropriate supercell which is 
depicted in Fig. [51 If we denote the primitive lattice 
vectors of graphene as a' and b', depicted as blue arrows 
in Fig. [2j the supercell vectors are given by a = 6a' and 
b = 2b' — a' (red arrows) thereby spanning a rectangular 
supercell containing 24 carbon atoms in the graphene 
layer. When the polymer chain direction is chosen to 
be parallel to b, its repeat unit is commensurate with 
the graphene basis vectors. Since we intend to study the 
behavior of a single PPP chain on graphene our choice 
of a = 6a' leads to a polymer chain separation which 
is large enough to prevent interaction between periodic 
replica of PPP chains. 

Fig. [3] shows the band structure of uncovered graphene 
in the above mentioned supercell containing 24 carbon 
atoms. This serves as a consistency test of our compu- 
tational approach since for uncovered graphene bands of 
the primitive hexagonal cell (compare Fig. [TJ) are sim- 



FIG. 2. (a) Perspective side view of the structural model 
describing poly(para-phenylene) (red) adsorbed on graphene 
(black). The adsorption height d is indicated, (b) Top 
view showing the unit cell vectors of the primitive, hexagonal 
graphene cell (blue arrows) as well as the rectangular super- 
cell spanned by vectors a and b (red arrows), (c) Brillouin 
zones corresponding to hexagonal graphene (black) as well as 
to the supercell (red). High symmetry points are indicated. 




FIG. 3. DFT-LDA (left) and GoWq (right) band structure 
of uncovered graphene calculated in the supercell (see text). 
The Fermi energy is set to zero. 
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FIG. 4. DFT-LDA (left) and Go Wo (right) band structure of 
poly(para-phenylene) (PPP). The Fermi energy is set to zero 
at the mid-gap energy. The frontier tt orbitals are indicated 
as explained in the text. 



ply folded into the Brillouin zone corresponding to the 
supercell. As shown in Fig. [5J:, the Dirac point, K, is 
folded to r, such that in the supercell there is a four-fold 
degeneracy at F leading to 4 touching cones. In T-X di- 
rection, just below the Fermi energy, there are two bands 
with slightly different slopes which are arising from the 
original Dirac cone along the directions KF and KM. 



B. Isolated PPP 

Before we investigate the adsorption of PPP on 
graphene, we present the electronic structure for an 
isolated PPP chain computed in the supercell intro- 
duced previously. The QP band structure obtained by 
G W @LDA is compared to LDA results in Fig. g] The 
valence band structure of PPP close to the Fermi level 
is characterized by two tt bands. The inter-ring bonding 
band, tti , is strongly dispersing from F down to Y with an 
LDA (G W ) band width of 3.54 eV (3.70 eV), while the 
inter-ring non-bonding band, 7T2, exhibits a much smaller 
dispersion along FY. Clearly, all bands along FX direc- 
tion, i.e. perpendicular to the polymer chain, exhibit a 
negligible dispersion reflecting the small inter-chain in- 
teraction. The frontier unoccupied bands 7r* and tt^ are 
the anti-bonding counterparts of m and 7T2, respectively 
The general trend of the self-energy corrections is to lower 
the occupied bands and raise the unoccupied bands rela- 
tive to the LDA results. This leads to an increase of the 
direct gap at the F point from 1.69 eV at the LDA level 
to 3.88 eV for the GqWq band structure. This finding is 
in accordance with previous calculations when con- 
sidering the fact that we neglect a possible torsion angle 




FIG. 5. Adsorption energy of PPP on graphene versus 
polymer-graphene distance d. Crosses are obtained by LDA, 
open triangles by using the PBE-GGA functional,— and filled 
circles by including an empirical van-der-Waals correction ac- 
cording to Grimme^ 



between adjacent phenyl rings and treat PPP as a per- 
fectly planar 7r-conjugated polymer. 



C. PPP adsorbed on graphene 

We now proceed to the main outcome of the paper, i.e. 
the electronic structure of PPP adsorbed on graphene. 
As a prerequisite, we first investigate the adsorption ge- 
ometry by varying the polymer-graphene distance and 
the adsorption site in DFT total-energy calculations. 
Fig. [5] shows the adsorption energy E a d as a function 
of adsorption height for the situation when the center of 
a PPP ring is positioned on top of a graphene carbon 
atom, as depicted in Fig. [2b, similar to the A-B stacking 
in bilayer graphene. The adsorption energy as well as 
the adsorption distance is very sensitive to the choice of 
the exchange-correlation potential. While LDA predicts 
an optimal adsorption distance of about 3.25 A with an 
adsorption energy of about 0.25 eV, GGAi2 results in 
almost no binding (0.02 eV) at the rather large distance 
of 4.5 A. This is, of course, indicative of a strong van- 
der-Waals contribution to the bonding^ For simplicity, 
we employ here an empirical scheme^ to correct for the 
missing dispersion forces in GGA, resulting in an adsorp- 
tion distance of d = 3.25 A and a binding energy of 
E a d = 0.39 eV as can be seen from Fig. [3] (red circles 
and line). When comparing the A-B type of adsorption 
site with an A-A type of adsorption position in which the 
hexagon of the polymer is on top of a graphene hexagon, 
we find an adsorption energy which is by 0.10 eV less fa- 
vorable. Based on the van-der-Waals corrected GGA^ 
the most favorable adsorption site of PPP on graphene is 
thus analogous to the A-B stacking in bilayer graphene 
and also in accordance with the adsorption position of a 
single benzene ring on graphite^ 
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FIG. 6. Plane-averaged total charge density and charge den- 
sity difference (top panel) as well as the electrostatic potential 
(lower panel) perpendicular to the graphene plane. The posi- 
tions of the graphene plane (diamond) and the polymer chain 
(triangle), the Fermi level, Ef, and the work function on bot- 
tom and top sides of the slab (f>t and (f>t are indicated. 

From the discussion above, it is evident that the bond- 
ing between the polymer and graphene is mainly due to 
van-der-Waals interactions. This is further emphasized 
by analyzing the charge-density rearrangements upon ad- 
sorption. To this end, we explore the charge density dif- 
ference defined in the following way. 

Ap = p - ( Ppo i+ p gr ). (3) 

Here p denotes the charge density of the combined sys- 
tem, while pp i and 

Pgr <irC the charge densities of iso- 
lated polymer and graphene, respectively. The plane- 
averaged charge density difference as well as the electro- 
static potential is depicted in Fig. [S] We can see that 
the interaction leads to regions of minor charge accumu- 
lation (A/? > 0) and charge depletion (Ap < 0) between 
the two constituents, i.e., to the formation of small sur- 
face dipoles. There is, however, no net charge transfer 
between the polymer and graphene. Consequently, the 
work-function modification due to the polymer adsorp- 
tion is negligible, thus the values on the bottom and top 
side of the slab, cf>b and <f>t respectively, agree with each 
other. These findings further prove that PPP is weakly 
physisorbed on graphene. 

In what follows, we present the band structure of the 
combined PPP-graphene system at the LDA as well as at 
the Go Wo level. In order to highlight the effect of image 
charges induced in the graphene layer on the electronic 
bands of PPP, we not only present data for the optimal 
adsorption height discussed above, but we also study the 
dependence of the band structure on the adsorption dis- 
tance d. We start the discussion by presenting the band 
structure of the adsorbate system at d = 4.0 A as shown 
in Fig. [7l Let us first focus on the LDA band structure 
depicted in the left panel of Fig. [7] Here, the solid black 
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FIG. 7. DFT-LDA (left) and Go Wo (right) band structure of 
PPP adsorbed on graphene for an adsorption height of 4 A 
shown as continuous black lines. The LDA band structure of 
an isolated PPP chain is shown as red, dotted line and used 
to identify PPP derived states in the combinbed system. The 
Go Wo band structure of the frontier PPP states is indicated 
by the red filled triangles, connecting lines serve as a guide for 
the eye. The GoWo enegies arising from graphene are shown 
as black filled circles. 



lines depict the bands of the combined system while the 
red, dashed lines indicate the bands of PPP obtained 
from a calculation for the isolated polymer. We observe 
that the LDA band structure of the combined system ap- 
pears as a mere superposition of the band strcutures of 
the isolated constituents, i.e. graphene and PPP. The 
relative alignment of the two subsystems' band struc- 
tures follows the simple rule of vacuum level alignment, 
also known as the Schottky-Mott limits Thus, also the 
HOMO-LUMO gap of the polymer, indicated by the red 
arrow at the T point, remains unaltered at the level of 
LDA-DFT. We note that at somewhat smaller polymcr- 
graphene distances, such as the equilibrium distance of 
3.25 A, 7T wave functions of the polymer and graphene 
start to overlap slightly leading to a weak hybridization 
between the polymer and graphene states. Correspond- 
ingly, also small charge re-arrangements and modifica- 
tions of the electrostatic potential due to the equilibra- 
tion of the chemical potentials of the two materials in 
contact take place as indicated in Fig. |5l However, these 
changes are of the order of 0.1 cV only. For weakly in- 
teracting systems, that are van-der-Waals bonded, we 
can summarize that the LDA band structure of the com- 
bined system is given by a superposition of the individual 
levels of the isolated subsystems, the polymer and the 
graphene sheet in our case, respectively. This is also the 
reason, why an almost substrate-independent behavior of 
the electronic properties of physisorbed organic molecules 
has been observed at the LDA-DFT level. 4 



Our findings for the the GqWq band structure, as de- 
picted in the right panel of Fig. are in stark contrast to 
this substrate- and adsorption-distance-independent be- 
havior of the LDA-DFT band structure. Here, the Go Wo 
band structure of the combined system is shown as black, 
filled circles where the frontier PPP bands are highlighted 
by the red triangles and continuous lines. When com- 
pared to the Go Wo results for the isolated PPP cahin 
(Fig. i}, a reduction of the HOMO-LUMO gap of PPP 
is evident. While the isolated chain exhibits a gap of 
3.95 eV, the value of the adsorbed PPP at a distance of 
4.0 A is reduced to 2.70 eV. This is a consequence of long- 
range correlation effects as a response of surface electrons 
to an added electron or hole in the polymer. This phe- 
nomenon is captured within the Go Wo approach by the 
screened Coulomb potential Wq. For molecules on sur- 
faces, this effect is particularly important as it contains 
the attractive interaction between the added electron or 
hole and its induced image charge. It can lead to a con- 
siderable reduction in the adsorbate's energy gap - about 
1.2 eV for the situation described above, where the band 
gap renormalization depends (i) on the polarizability of 
the surface and (ii) on the adsorbate's distance from the 
surface. Both effects arc absent in the DFT electronic 
structure at the LDA or GGA level due to the locality 
of exchange-correlation potential, and this deficiency can 
also not be cured by using hybrid functionals4. 
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FIG. 8. The energy of the HOMO (filled triangles) and LUMO 
(open squares) of PPP as a function of the adsorption distance 
d on a reciprocal scale. The left panel shows LDA results 
while the right panel displays corresponding GoWo values. 
The lines represent various electro-static models as described 
in the text. Note that energies are measured with respect to 
the vacuum level. 



D. Adsorption distance dependence 

To further emphasize the physical origin of the band- 
gap reduction upon adsorption, we study its dependence 
on the adsorption distance d. The results are depicted 
in Fig. [8] in which we plot the energies of the HOMO 
and LUMO of PPP at the V point as a function of PPP- 
graphene distance d, where the left (right) panel displays 
LDA (Go Wo) results. Note that the limiting values for 
d = oo are taken from calculations for an isolated PPP 
chain. Let us first discuss the distance dependence for 
moderately large values of d > 5 A. Here, the HOMO and 
LUMO energies, i.e., the ionization and electron affini- 
ties levels of PPP, are independent of d in LDA while 
they show strong and opposite trends in the Go Wo cal- 
culations leading to the above mentioned reduction of 
the band gap. This asymptotic distance-dependence can 
be understood by employing a classical image-potential 
model of the form22 



Here, q' = — q(e — l)/(e+ 1) is the induced image charge, 
where e is the relative dielectric constant of the polar- 
izable medium, and do denotes the effective position of 
the image plane. It describes the electrostatic interac- 
tion between charges (added electrons or holes) above 
the surface and the polarization charge below the sur- 
face. This interaction also constitutes the basic term 



of the Go W) self-energy^i and is consistent with the 
physical picture outlined above, i.e., that the level re- 
normalizations are determined by long-range Coulomb 
interactions which decay very slowly with distance. The 
solid, black lines in Fig. [5] display the expected image- 
potential corrections according to Eq. ((4]) when inserting 
an effective image plane position do = 1.0 A and a dielec- 
tric constant of e = 3. The former value is in accordance 
with the image plane position computed for benzene ad- 
sorbed on graphite,— while the latter agrees with a re- 
cent result by Wchling et al. who calculated the dielec- 
tric constant of graphene to be 2.4.— Using these model 
parameters, the image model (solid line) only roughly 
follows the Go W) values (symbols). It underestimates 
the level re-normalizations for large distances, particu- 
larly so for the HOMO, and does not capture the fact 
that the GqWq electron affinity level seems to saturate at 
polymer-graphene distances smaller than 4 A, and also 
can not explain the more complicated behavior of the 
ionization potential at small d. 

We further investigate possible reasons for this devia- 
tion. First, a better agreement between the image model 
(|4|) and Go Wo results for large distances can be achieved 
when allowing the position of the effective image plane 
to move further away from the graphene plane, e.g., to 
do = 1.5 A and also assuming a slightly bigger e of 4 
(black, dashed lines). These numbers lie still within the 
range of values reported for similar surfaces in an ear- 
lier work^ and improve the agreement for separations 
d > 5 A, while they clearly worsen the description for 
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smaller distances. It is evident that adjusting do and 
e will not lead to a satisfactory description of the data 
points over the full range of distances due to the leveling 
off at small polymer-graphene separations. This indicates 
another physical mechanism that becomes important for 
smaller distances. It has been noted earlier,— that the 
simple model given in Eq. Q based on a point-like charge 
does not take into account the polarizability of the ad- 
sorbate, an effect which should be more pronounced at 
small distances. To estimate the magnitude and direction 
of such a polarization effect, we have extended the stan- 
dard image potential model as to to include the polar- 
izability of the adsorbate system in a simplified manner 
(see Appendix iBl for details). To this end, we consider 
that the image charge q' induces a dipole moment in the 
adsorbate system which in turn gives rise to a dipole 
field acting on the physical charge. This model (orange 
line) qualitatively follows the computed Go Wo values. It 
agrees with the standard model at large d but deviates 
from it at small d due to a correction of opposite sign 
whose leading l/<i 3 -term is proportional to the polariz- 
ability of the adsorbate. Thus, the level-renormalizations 
level off at small d in accordance with the Go Wo. While 
our extended model provides physical insight as to why 
the standard electro-static model starts to break down, 
we emphasize that it is still too crude to capture details 
such as the shape of the adsorbate or the anisotropy of 
its polarizability. Thus, it may serve as a tool for qualita- 
tive understanding but not for a quantitative predictions. 
Moreover, there are also other effects which may cause 
the standard image charge model of Eq. (TJ| to fail for 
small adsorption distances. For instance, for polymer- 
graphene distances smaller than 4.5 A, the added electron 
(hole) into the 7r* (m) state at T starts to extend beyond 
the position of the fictitious image plane which may also 
lead to a deviation from the 1/d behavior. Such small 
hybridization effects between polymer and graphene are 
exemplified by the minor, but clearly visible, dependence 
of the HOMO and LUMO LDA-energies (left panel of 
Fig. HJ. Finally, to obtain the full electronic band struc- 
ture, one needs to determine the individual self-energies 
for each k point and band, which goes far beyond the 
capability of such model. 
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FIG. 9. Comparison of the Go Wo gap of graphene at the T 
point using the conventional sum-over-empty-states approach 
(circles and solid black line) and the energy effective technique 
(red line). 

proximation, we obtain a considerable reduction of its 
HOMO-LUMO gap upon adsorption even for large dis- 
tances from graphene where the wavefunction overlap be- 
tween graphene and PPP is negligible. We find that a 
classical image-potential model in its standard form de- 
scribes the G Wo HOMO and LUMO energies of PPP 
only for fairly large distances, while we observe some de- 
viations from the expected 1/d dependence for polymer- 
graphene separations smaller than 4.5 A down to the 
equilibrium van-der-Waals bond distance of 3.25 A. By 
incorporating the polarizability of the adsorbate into the 
model, we are able to qualitatively improve it also for 
shorter polymer-graphene distances, thereby identifying 
the most-likely physical mechanisms for the substrate- 
induced level renormalizations of the polymer HOMO 
and LUMO close to equilibrium bonding distances. How- 
ever, for quantitative description Go Wo calculations are 
necessary to properly predict the electronic structure of 
the adsorbates close to surfaces. 
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IV. SUMMARY 
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structure calculations. Analysis of the charge density 
shows the formation of weak interface dipoles, but no 
net charge transfer between PPP and graphene. At 
the level of density-functional theory within the local- 
density approximation, we find that the adsorption of 
PPP on graphene docs not alter the PPP band struc- 
ture compared to an isolated PPP chain. However, by 
incorporating many-body effects within the Go Wo ap- 
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Appendix A: Convergence tests for monolayer 
graphene 

As discussed in Section [Til Go Wo calculations typically 
involve time-consuming summations over empty states. 
Here, we demonstrate for the case of graphene the con- 
vergence of the QP energies with respect to the number 
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FIG. 10. The convergence of the Go Wo band gap at F point 
for monolayer graphene when using a cut-off version of the 
Coulomb potential (filled symbols) compared to the plain 
Coulomb potential (open symbols). The left panel shows the 
convergence as a function of vacuum-layer thickness c for a 
fixed fc-point sampling of 30 x 30 x 1 while the right panel dis- 
plays the convergence as a function of fc-points {Nk m X Nk x x 1 
meshes) for a fixed vacuum layer thickness of c = 15.9 A. 



of empty states and compare with results from the energy 
effective technique in which all summations are restricted 
to occupied states only. This is shown in Fig.|H]for the QP 
energy differences between the highest occupied and the 
lowest unoccupied band at the T point. Note that for this 
convergence test, we have used a moderately dense sam- 
pling of only 6x6x1 fc-points. One can see that the QP 
energy differences obtained by both methods converge 
to the same value. In view of this very good agreement 
and the computational efficiency of the EET, we apply 
this technique in all GqWq calculations presented in this 
paper. 

As outlined in Sec. HH the nonlocal nature of the self- 
energy leads to a slow convergence of results as a func- 
tion vacuum layer size in slab geometries which can be 
mediated by using a truncated Coulomb potential. The 
left panel of Fig. [TU] shows the convergence of the Go Wo 
gap at the T point of hexagonal graphene. The open 
squares display results for the unmodified Coulomb po- 
tential while the filled squares are values obatincd by 
using a cut-off Coulomb potential for slab geometries 
according to Ismail-BeigiJ^ Clearly, the quasi-particle 
gap is converged at a graphene layer separation of about 
c = 20 A when using the modified Coulomb potential, 
while it approaches a converged value only via a 1/c de- 
pendence in the unaltered Coulomb potential (note the 
reciprocal abscissa in Fig. [TO]). For instance, at 20 A, 
the error would amount to about 0.2 eV this particu- 
lar quasi-particle energy difference. It should be noted, 
however, that the fc-point convergence worsens when uti- 
lizing a cut-off Coulomb potential. This is visualized in 
the right panel of Fig. [TU1 where we plot the Go Wo gap 
of graphene at T as a function of Nk x which defines the fc 




FIG. 11. Schematics showing the image-potential model. A 
charge q at position z above a dielectric with dielectric con- 
stant e filling the half space z < induces an image charge 
q' = — which, in turn, induces a dipole with dipole mo- 
ment p z at the location d of the adsorbate which is modeled 
by a homogeneous polarizable sphere of radius R. The electric 
fields arising from the image charge (red line) as well as from 
the image-charge-induced dipole (orange line) are illustrated 
for d — 4 A, a = 15 a.u. and R = 2.5 a.u. 



mesh (Nk x x -^fet x 1 points). From the analysis of Fig.fTUl 
we conclude that a vacuum layer thickness of 20 A and 
a fc-mesh of 36 x 36 x 1 should be sufficient to converge 
quasi-particle energies to within 0.1 eV when employing a 
cut-off Coulomb potential. For the supercell used in the 
calculation of the polymer and the combined polymer- 
graphene system, this fc-mesh translates into a mesh of 
5 x 17 x 1 which we have utilized throughout the paper 
for all computations involving the supercell. 



Appendix B: Image-potential model 

The classical image-potential model leading to Eq. (HJ) 
assumes a point charge in front of a semi-infinite (z < 0), 
polarizable medium described by the dielectric constant 
e. A charge q which resides a distance z = d above the 
dielectric polarizes the dielectric, an effect which is ac- 
counted for by an image charge q' = —q(e — l)/(e + 1) 
appearing at z = —d. When integrating the attractive 
force between the charge and its image F(d) = jjjyz from 
d to infinity, the energy correction given in Eq. (j4]) is ob- 
tained. Note that we use atomic units throughout, i.e., 
we set e 2 /47reo — > 1. 

Let us now extend this simple point-charge description 
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by allowing the adsorbatc system to be polarizable. As 
indicated in Fig. [TTJ we imagine an adsorbed atom at a 
fixed distance d above the semi-infinite dielectric. When 
adding a charge q at distance z, the induced charge q' at 
— z will now polarize the atom thereby inducing a dipole 
moment p z of strength 



Pr, 



aE z 



(Bl) 



Here, a denotes the polarizability of the adsorbate, and 
E z is the z-component of the electric field at the location 
of the physical charge due to the image charge q' , thus 



E z = 



{z + d) 2 ' 



(B2) 



This induced dipole gives rise to an electric field at z 
which exerts a force on the charge q in addition of the 
image force due to q' . Combining the above expressions 
one finds for the z-component of the total force 



F z = qq' 



1 



2a 



(2z) 2 (z + df\z - d\ ? - 



(B3) 



Here, the first term (red line in Fig. ITTj) describes the 
attraction of the charge q by the image charge q' , while 
the second term (orange line), which is proportional to 
the adsorbate's polarizability a, arises due to the image- 
chargc-induccd dipole field. When integrating the force 
from d to oo, the first term gives Eq. @ while the sec- 
ond term produces an expression whose leading term is 
proportional to 1/d 3 . Two things should be noted at this 
stage. First, the induced dipole at z = d in turn results 
in an image dipole at z = —d which, when taken into 



account in self-consistent manner, exerts a force on the 
charge q at z. However, this effect which is proportional 
to a 2 modifies the final result only at leading order of 
1/d 6 , i.e., at very short distances due to the faster de- 
cay of the dipole compared to the monopole field. As 
a second note, the dipole field proportional to l/\z — d\ 3 
would lead to a diverging force if the charge q approaches 
the position of the atom z = d. In order to circumvent 
this problem arising from a too simple model, we rep- 
resent the polarizable adsorbate by a homogeneous, di- 
electric sphere of radius R exhibiting the polarizability 
a. Then, the field inside the sphere is behaves regu- 
larly at z — > d and the force may be readily integrated 
yielding the energy rcnormalization Vi m as a function of 
adsorbatc-substrate separation d 



V im (d) = qq' 



1 3a 1 

Id~8Rd 3 + { dA> 



(B4) 



The first term gives again the standard point-chargc- 
result, while the remaining terms starting with 1/d 3 are 
due to polarization effects of the adsorbate and conse- 
quently contain the polarizability a and the size R of the 
adsorbate system. When further allowing for a shift of 
the image plane to the position do, values of a = 15 a.u. 
and R = 2.5 a.u. result in the orange curve shown in the 
right panel of Fig. [5] We note that the atomic polar- 
izability of an isolated carbon atom is about 11 a.u.^2 
and emphasize that these numerical values merely serve 
to underpin the main physical effect rather than to cre- 
ate quantitative predictions which would require a more 
refined electro-static model which takes into account the 
actual shape, the inhomogeneity as well as the anisotropy 
of the adsorbate system. 
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